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ABSTRACT 


We  have  used  the  Navy  Prototype  Optical  Interferometer  (NPOI)  to  obtain  the  first  multichannel  optical 
aperture  synthesis  images  of  a  star.  We  observed  the  spectroscopic  binary  i Ursae  Majoris  at  6  to  10 
milliarcseconds  separation  during  seven  nights,  using  three  interferometric  baselines  and  19  spectral 
channels  (XX  520—850  nm)  of  the  NPOI.  After  editing,  a  typical  90  sec  scan  yielded  fringe  visibilities  at 
50  spatial  frequencies  and  closure  phases  at  15  wavelengths.  Three  to  five  scans  were  obtained  each  night. 
The  separations  and  position  angles  are  in  good  agreement  with  the  visual  orbit  obtained  with  the  Mark  III 
interferometer  (Hummel  et  al.  1995,  AJ,  110,  376)  but  show  small  systematic  difference  that  can  be  used 
to  improve  the  orbit.  The  closure  phase  data  provide  a  sensitive  measure  of  the  magnitude  difference 
between  the  components.  These  results  demonstrate  the  power  of  broad-band  interferometric  observations 
for  fast  imaging  and  the  utility  of  vacuum  delay  fines  for  simultaneous  observations  over  a  wide  band.  These 
observations  are  the  first  to  produce  simultaneous  visibilities  and  closure  phases  with  a  separate-aperture 
optical  interferometer,  and  the  second  to  produce  closure  phase  images,  following  the  results  from  COAST 
reported  by  Baldwin  et  al.  (1996,  A&A,  306,  L13).  The  angular  resolution  here  is  the  highest  ever  achieved 
at  visual  wavelengths,  exceeding  by  an  order  of  magnitude  the  best  thus  far  achieved  by  any  single-aperture 
optical  telescope.  We  generated  complex  visibilities  and  closure  phases  (the  data  types  commonly  used  in 
radio  interferometry)  from  the  optical  data  and  used  standard  radio  interferometry  techniques  to  produce 
these  images.  However,  the  fundamental  observables  of  optical  interferometry,  the  squared  visibility 
amplitude  and  the  closure  phase,  require  the  development  of  new  analysis  techniques.  ©  1997  American 
Astronomical  Society.  [S0004-6256(97)01709-3] 


1.  INTRODUCTION 

Atmospheric  turbulence  limited  the  resolution  of  tele¬ 
scopes  at  visual  wavelengths  until  early  in  the  twentieth  cen¬ 
tury,  when  A.  A.  Michelson  demonstrated  that  interferom¬ 
etry  had  the  potential  of  measuring  the  exceedingly  small 
angular  diameters  of  stars  (Michelson  &  Pease  1921).  Mich¬ 
elson  also  realized  that  measurements  of  fringe  visibilities 
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could  be  used  to  image  stellar  surfaces.  The  use  of  interfer¬ 
ometry  to  measure  orbits  of  binary  stars  also  began  early  in 
this  century  (Anderson  1920).  Brown  et  al.  (1974)  used  in¬ 
tensity  interferometry  to  measure  stellar  angular  diameters. 
Speckle  interferometry,  a  somewhat  different  technique,  has 
also  been  used  to  measure  binary  orbits. 

Atmospheric  turbulence  limits  interferometry  as  well,  by 
rapidly  shifting  the  interference  fringes.  After  Johnson  et  al. 
(1974)  and  Labeyrie  (1975)  showed  that  the  fringe  amplitude 
can  be  measured  with  an  interferometer  consisting  of  me¬ 
chanically  separate  telescopes,  Shao  &  Staelin  (1977)  dem¬ 
onstrated  that,  with  computer-controlled  delay  lines  and 
high-speed  detectors,  the  optical  fringe  phase  can  be  stabi- 
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lized  and  tracked  for  hours.  The  Mark  III  optical  interferom¬ 
eter  on  Mt.  Wilson  (Shao  et  al.  1988)  advanced  this  technol¬ 
ogy  and  demonstrated  the  capabilities  of  a  two-element 
interferometer  for  measuring  stellar  diameters  (e.g.,  Hutter 
et  al.  1989;  Mozurkewich  et  al.  1991,  1997)  and  positions 
(e.g.,  Hummel  et  al.  1994),  as  well  as  binary  orbits  (e.g., 
Hummel  et  al.  1995). 

The  next  step  in  optical  interferometry  is  the  extension  to 
multiple  baselines,  to  which  the  phase-closure  techniques 
first  espoused  by  Jennison  (1958)  and  subsequently  devel¬ 
oped  for  very  long  baseline  radio  interferometry  (VLBI)  can 
be  applied.  In  general,  closure  phase  is  needed  to  obtain  un¬ 
ambiguous  images  from  an  interferometer.  A  two- 
dimensional  multibaseline  array  with  baseline  lengths  of  100 
m  or  more  could  image  stars  with  sub-milliarcsecond  reso¬ 
lution.  This  paper  describes  a  significant  milestone  achieved 
in  the  construction  of  such  an  array,  the  Navy  Prototype 
Optical  Interferometer  (NPOI),  a  joint  project  of  the  US  Na¬ 
val  Research  Laboratory  and  the  US  Naval  Observatory  in 
cooperation  with  Lowell  Observatory. 

We  report  the  first  images  obtained  with  the  NPOI,  made 
with  simultaneous  multichannel  closure  phase  and  visibility 
amplitude  data.  We  observed  the  spectroscopic  binary  £'  Ur- 
sae  Majoris  [Mizar  A,  mv=  2.27,  P  =  20'?54,  a  =  9.64  mil- 
liarcseconds  (mas)  (Hummel  et  al.  1995)]  on  seven  nights 
between  1996  May  1  and  June  4.  Processing  our  closure 
phases  and  visibility  amplitudes  with  standard  interferomet¬ 
ric  image  synthesis  techniques,  we  have  achieved  a  spatial 
resolution  of  «=3  mas.  Our  achieved  resolution  exceeds  by 
factors  of  two  to  three  the  best  resolution  theoretically  ob¬ 
tainable  by  the  current  generation  of  monolithic  optical  tele¬ 
scopes.  The  closure  phases  produce  a  sensitive  measure  of 
the  brightness  ratio  of  the  two  components.  We  show  here 
that,  by  taking  data  simultaneously  in  multiple  spectral  chan¬ 
nels  (in  these  images,  we  use  19  channels  spanning  a  factor 
of  1.6  in  spatial  frequency),  we  greatly  increase  the  speed 
with  which  images  can  be  made. 

The  observations  described  here  produced  only  the  sec¬ 
ond  closure  phase  measurements  ever  made  with  a  separate- 
element  optical  interferometer.  They  are  the  first  observa¬ 
tions  with  simultaneous  closure  phase  and  visibility 
amplitude  measurements,  and  the  first  with  multibaseline 
data  taken  in  more  than  one  spectral  channel.  The  first 
closure-phase  measurements,  in  observations  of  the  56  mas 
binary  a  Aurigae  (Capella),  were  reported  by  the  COAST 
group  at  Cambridge  University  (Baldwin  et  al.  1996). 

2.  THE  NPOI 

The  NPOI,  located  at  the  Lowell  Observatory  site  on 
Anderson  Mesa,  near  Flagstaff,  Arizona,  is  a  long-baseline 
optical  interferometer  which  will  consist  of  an  astrometric 
sub-array  (four  elements)  and  an  imaging  sub-array  (six  ele¬ 
ments)  when  construction  is  complete  (Armstrong  et  al. 
1997).  The  purpose  of  the  NPOI  is  to  determine  the  precise 
positions  of  bright  stars  and  to  image  the  surfaces  and  cir- 
cumstellar  environments  of  nearby  stars  and  stellar  systems. 
The  array  elements  are  50  cm  siderostats;  however,  only  12.5 
cm  of  the  aperture  is  currently  used,  a  limit  set  by  the  size  of 


the  feed  system  optics.  Beam-compressing  telescopes,  which 
will  increase  the  useable  apertures  to  35  cm,  are  currently 
being  built  for  the  astrometric  sub-array. 

The  NPOI  uses  32  spectral  channels  covering  a  wide  band 
(XX  450—  850  nm)  to  increase  its  sensitivity  and  its  range  of 
spatial  frequencies,  and  to  help  determine  the  effect  of  atmo¬ 
spheric  longitudinal  dispersion  on  the  fringe  positions.  Each 
channel  corresponds  to  a  spatial  frequency  of  D/X,  where  D 
is  the  baseline  length  and  X  is  the  mean  wavelength.  In  ad¬ 
dition,  the  NPOI  uses  vacuum  delay  lines  to  ensure  that 
fringes  over  the  entire  band  can  be  observed  simultaneously. 
Without  vacuum  delay  lines,  the  large  differences  in  air  path 
length  between  the  individual  light  paths  and  the  variation  in 
the  index  of  refraction  of  air  over  the  visual  range  produce  a 
large  variation  of  fringe  phase  with  wavelength,  precluding 
broadband  observations  that  extend  into  the  blue  part  of  the 
visual  spectrum.  The  instrument  incorporates  active  group- 
delay  fringe  tracking  and  rapid  tip-tilt  star  tracking.  Ava¬ 
lanche  photodiodes  are  used  in  both  the  fringe  detection  and 
star  tracking  systems. 

The  first  single-baseline  stellar  fringe  measurements  at  the 
NPOI  were  obtained  in  1994  October.  After  numerous  equip¬ 
ment,  infrastructure,  and  software  upgrades,  including  the 
addition  of  the  third  siderostat,  the  first  three-baseline  stellar 
fringe  observations,  including  closure-phase  measurements, 
were  made  on  1996  March  18. 

As  with  the  Mark  III  interferometer,  the  NPOI  modulates 
the  delay  on  each  baseline  by  a  small  number  of  wavelengths 
and,  for  each  spectral  channel  and  baseline,  synchronously 
measures  the  photon  count  rates  in  n  bins  per  wavelength  of 
modulation.  (Currently,  n  =  8  for  the  NPOI;  for  the  Mark  III, 
n  =  4.)  We  generate  the  complex  Fourier  transform  of  this  set 
of  bins  and  select  the  component  corresponding  to  a  fre¬ 
quency  of  one  fringe  per  n  bins.  The  real  and  imaginary  parts 
of  this  component  are  X  and  Y .  respectively.  The  unbiased 
estimator  for  the  squared  visibility  amplitude  V 1  for  one 
baseline  is  given  by 


4 

(X2+Y2-v2n) 

sinc2(l/n) 

(N)2 

Here,  N  is  the  total  number  of  photons  in  a  frame  of  n  bins, 
and  sinc(x)  =  sin(ire)/77x.  Subtracting  <j2n,  the  variance  of  N, 
from  X2+  Y2  collects  for  the  bias  due  to  noise  (see  Tango  & 
Twiss  1980,  who  treat  the  case  of  shot  noise,  in  which 
crN=N).  We  calculate  V2  rather  than  |L|  because  the  bias 
correction  for  |F|  depends  on  the  value  of  j  V\  itself.  The 
numerator  and  denominator  are  computed  separately  for  ev¬ 
ery  2  ms  sample  interval  and  averaged  for  1  s.  The  1  s 
records  are  edited  and  further  averaged  over  the  length  of  the 
scan.  The  complex  Fourier  components  X+iY  from  the 
three  baselines  are  multiplied  together  every  2  ms  to  give  the 
complex  triple  product.  We  average  the  triple  product  coher¬ 
ently,  since  the  atmospheric  phase  noise  cancels.  The  phase 
of  the  triple  product  is  the  closure  phase,  <pcl . 

The  calibration  of  the  observed  average  V2  and  <pA  is 
based  on  unresolved  single  stars.  For  unresolved  sources,  the 
theoretical  V2  and  4>o\  are  one  and  zero,  respectively.  The 
deviations  of  the  measured  V2  and  tf>A  of  the  unresolved 
calibrators  from  these  values  can  be  written  as  low-order 


©  American  Astronomical  Society  •  Provided  by  the  NASA  Astrophysics  Data  System 


1997AJ. 


1223  BENSON  ET  AL.:  f1  UMa 


1223 


_Q 

‘(f) 

> 

~o 

<D 

L_ 

o 

=3 

CT 

00 


1 .0 

i 

- ■ - , - 1 - ,  t  i  |  i  > 

+ 

0.8 

+ 

+ 

+ 

0.6 

+ 

+ 

+ 

+ 

0.4 

+ 

+  + 

0.2 

" 

+ 

+  +  + 

++ 

: 

0.0 

-  ( 

_ . _ . _ 1 _ . - . - . - 1 - 1 - . — 

.  1 

400  600  800  1000 

Wavelength  [nm] 


Fio.  1.  Squared  visibility  amplitude,  V2,  and  closure  phase,  <f>d ,  vs  wavelength  for  the  first  scan  of  1996  May  1.  The  V2  shown  are  from  the  spectrometer  for 
the  EW  baseline;  data  from  bad  detectors  and  from  the  channel  contaminated  by  scattered  HeNe  laser  metrology  light  were  deleted.  The  closure  phase 
combines  data  from  all  three  spectrometers,  one  per  baseline.  A  bad  channel  in  any  spectrometer  precludes  measuring  <f>d  at  that  wavelength;  thus  more 
channels  are  missing  from  these  data  than  from  the  EW  V2  data.  For  both  V2  and  <f>c] ,  wavelengths  shortward  of  X520  nm  were  deleted  because  of  low 
signal-to-noise  ratio.  The  formal  uncertainties  in  V1  are  smaller  than  the  symbols,  but  the  calibration  uncertainties  are  ~  ±  10%.  For  d>(] ,  the  formal 
uncertainties  are  again  smaller  than  the  symbols,  but  the  calibration  uncertainties  are  ~±2°. 


polynomial  functions  of  atmospheric  and  instrumental  vari¬ 
ables.  The  polynomial  coefficients  are  determined  from  the 
calibrator  scans  and  subsequently  applied  to  the  data  of  both 
the  program  sources  and  the  calibrators.  The  deviation  be¬ 
tween  the  calibrated  data  of  the  calibrator  stars  and  their 
expected  values  is  a  measure  of  the  accuracy  of  our  calibra¬ 
tion  procedure. 


3.  THE  DATA  AND  RESULTING  IMAGES 

We  used  the  center  (C),  east  (E),  and  west  (W)  siderostats 
of  the  astrometric  sub-array  to  observe  £*  Ursae  Majoris  on 
1996  May  1-4  and  29,  and  June  1  and  4.  Several  calibration 
stars  were  observed  each  night,  including  a  nearby  calibrator, 
y  UMa.  The  baseline  lengths  were  18.9,  22.2,  and  37.5  m  at 
azimuths  —67!  5,  63!  6,  and  86!  0,  respectively.  Between 
three  and  five  scans  (on-source  integrations,  typically  90  s 
each)  were  recorded  on  '(l  UMa  each  night,  with  a  similar 
number  of  scans  on  y  UMa.  We  deleted  12  channels  at  the 
blue  end  (A <520  nm)  due  to  low  signal-to-noise  ratio.  We 
also  removed  a  small  number  of  channels  with  bad  detectors, 
as  well  as  the  channel  contaminated  by  the  HeNe  delay-line 
metrology  laser.  Typically,  each  scan  produced  measure¬ 
ments  of  V2  at  *=50  spatial  frequencies.  Measuring  <f>cl  at  a 
given  wavelength  requires  valid  data  in  the  corresponding 
channel  from  all  three  baselines.  With  this  requirement,  the 
typical  scan  produced  15  measurements  of  cf>d. 

For  the  images  presented  here,  a  single  calibration  factor 
for  V2  was  derived  from  the  calibrator  data  for  each  channel, 
baseline,  and  night  independently.  For  (f>d ,  a  low-order  poly¬ 
nomial  versus  time  was  fitted.  These  calibrations  used  all  the 
calibrators  observed  during  a  given  night.  We  are  developing 
more  sophisticated  calibration  methods  based  on  procedures 
used  for  the  Mark  III  interferometer  (Shao  et  al.  1988;  Mo- 
zurkewich  et  al.  1991;  Armstrong  et  al.  1992). 


The  uncertainties  in  the  resulting  V2  and  (f>c]  are  domi¬ 
nated  by  calibration  uncertainties;  the  formal  errors  are 
~±2%  for  V2  and  ~±0!5  for  (f>d,  but  the  calibration  un¬ 
certainties  are  ~±10%  for  V2  and  —  ±2°  for  cf>d  ■  The 
greatest  effect  of  these  uncertainties  is  on  the  magnitude  dif¬ 
ference  Am;  the  separations  and  position  angles  of  the  binary 
components  are  determined  primarily  by  the  location  of  the 
maxima  and  minima  of  V2  in  the  (u,v)  plane,  which  are  only 
weakly  dependent  on  the  calibration. 

Figure  1  shows  the  E-W  V 2  data  and  the  <f>d  data  for  one 
of  the  May  1  scans.  The  sinusoidal  variation  of  V2  with 
increasing  wavelength  (i.e.,  decreasing  (u,v)  spacing)  shows 
that  the  binary  is  clearly  resolved  on  this  baseline.  The 
change  in  <t>c\  of  nearly  180°  at  the  wavelength  of  the  mini¬ 
mum  in  V2  indicates  that  the  two  components  of  the  binary 
are  of  nearly  equal  brightness. 

Because  the  components  of  £*  UMa  are  very  similar  (Pet¬ 
rie  1939;  Fehrenbach  &  Prevot  1961;  Hummel  et  al.  1995), 
the  image  structure  is  wavelength  independent,  and  one  can 
treat  the  spectral  channels  as  separate  (u,v)  samples  of  a 
single  image.  The  quality  of  such  a  multichannel  optical  ap¬ 
erture  synthesis  image  is  much  better  than  that  of  a  single¬ 
channel  image,  since  the  factor  of  —1.6  in  wave  number 
range  across  our  observed  band  expands  the  spatial  fre¬ 
quency  coverage  by  the  same  factor.  The  orbital  period  is 
20'? 54,  so  the  orbital  motion  is  small  during  one  night.  For 
that  reason,  we  treat  all  of  a  given  night’s  scans  together, 
which  makes  an  additional  significant  improvement  in  («,i>) 
coverage.  Figure  2  shows  the  (u,v)  coverage  for  the  four 
scans  taken  on  May  1  treated  as  a  single  (u,v)  data  set. 

We  used  the  Astronomical  Image  Processing  System 
(AIPS)  to  make  one  image  from  the  data  for  each  night.  For 
this  purpose,  we  constructed  complex  visibilities  for  each 
baseline  and  spectral  channel  by  assigning  (V2)1/2  to  the  vis¬ 
ibility  amplitudes  and  zero  to  two  of  the  visibility  phases. 
We  then  assigned  the  calibrated  4>a  to  the  phase  of  the  vis- 
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Fig.  2.  Spatial  frequency  coverage  for  the  binary  Mizar  A  on  1996  May  1 
with  the  center,  east,  and  west  (“C,”  “E,”  “W”)  elements  of  the  astro- 
metric  sub-array  of  the  NPOI.  Each  baseline  produces  data  from  32  spectral 
channels  and  thus  results  in  a  ray  in  the  (u,v)  diagram.  Four  scans  of  (u.v) 
data  are  shown  here,  so  each  baseline  produces  four  rays.  The  longest  (u,v) 
spacings  are  from  the  east-west  (“EW”)  baseline.  “CW”  (upper  four  rays) 
and  “CE”  (lower  four  rays)  indicate  data  from  the  two  shorter  baselines. 


ibility  on  the  E-W  baseline.  This  procedure  produced  a  set  of 
125  to  180  visibilities  (depending  on  the  number  of  scans 
and  of  available  t/>cl  measurements)  from  which  to  generate 
an  image  for  each  night.  Standard  phase  self-calibration 
methods  employed  in  AIPS  readjust  the  baseline  phases 
while  keeping  <pcl  constant.  We  processed  the  data  with  uni¬ 
form  weighting  and  a  cell  size  of  1  mas,  and  produced 
256X256  pixel  maps.  For  the  May  1  data,  the  resulting  syn¬ 
thesized  beam  was  5.3X  1.7  mas  at  position  angle  6°.  This 
synthesized  beam  size  was  typical  for  all  seven  nights  of 
observation. 

We  CLEANed  each  image,  used  the  CLEAN  components 
as  a  model  for  self-calibrating  the  phases,  and  re-imaged. 
The  resulting  image  for  the  May  1  data,  after  a  few  itera¬ 
tions,  is  shown  in  Fig.  3.  The  dynamic  range  (peak  to  rms)  of 
this  image  is  >100:1.  Multichannel  optical  aperture  synthe¬ 
sis  images  for  the  other  six  nights,  produced  in  a  similar 
manner  to  the  May  1  image,  are  shown  in  Fig.  4.  The  binary 
is  clearly  resolved  in  all  seven  images,  and  the  relative  mo¬ 
tion  of  the  components  is  evident. 

We  fitted  Gaussian  intensity  profiles  to  the  images  to  de¬ 
termine  their  relative  positions.  Figure  5  shows  these  relative 
positions  from  the  seven  nights  of  observation,  as  well  as  the 
visual  orbit  (semimajor  axis  a  =  9.6  mas,  period  P  =  20? 54) 
determined  with  the  Mark  III  optical  interferometer  (Hum¬ 
mel  et  al.  1995).  The  measured  NPOI  positions  are  con¬ 
nected  by  lines  to  the  positions  predicted  from  the  Mark  III 
orbit.  The  uncertainty  ellipses  are  one-fifth  the  FWHM  of  the 
synthesized  beam,  a  standard  criterion  used  in  radio  VLBI 
imaging.  The  agreement  with  the  Mark  III  orbit  is  good,  but 
shows  slight  systematic  differences,  which  may  result  from  a 
gap  in  orbital  phase  coverage  in  the  Mark  III  data. 

We  also  note  that  the  shape  of  the  <£cl(X)  curve  and  the 
size  of  the  jump  in  4>a  are  sensitive  indicators  of  Am.  For 
Am  =  0,  <pc\(X)  is  a  step  function  of  amplitude  180°,  with 


Fig.  3.  Image  of  Mizar  A  from  1996  May  1.  The  coordinates  are  offsets  in 
milliarcseconds  between  the  components.  Contour  levels  are  —  1  %,  —0.5%, 
0.5%,  1%,  2%,  5%,  10%,  20%,  50%,  and  80%  of  the  peak  flux.  Both  stars 
have  angular  diameters  ^0.8  mas,  and  thus  are  unresolved  in  these  images. 
The  restoring  beam  was  a  circular  Gaussian  of  3  mas  diameter  (FWHM). 

the  step  at  the  wavelength  of  the  minimum  of  the  triple  am¬ 
plitude.  (For  the  scan  shown  in  Fig.  1,  this  occurs  at  the 
minimum  of  V 2  for  the  EW  baseline.)  For  Am  ¥=  0,  the  size  of 
the  jump  is  less  than  180°,  and  the  comers  of  the  step  func¬ 
tion  become  rounded.  For  these  data,  we  have  determined  a 
preliminary  value  for  Am  of  OfOlLOfOl,  in  agreement 
with  earlier  estimates  that  Am  is  close  to  zero  (Petrie  1939; 
Fehrenbach  &  Prevot  1961;  Hummel  et  al.  1995).  Our  mea¬ 
surement  of  Am  applies  near  X700  nm,  where  we  see  the 
jump  in  (f)cl .  Since  the  two  components  of  'C  UMa  are  very 
similar,  Am  should  be  close  to  zero  across  the  entire  band, 
but  measuring  Am  as  a  function  of  X  requires  more  data. 

The  images  generated  with  ABPS,  however,  show  Am  val¬ 
ues  of  up  to  O’M;  only  the  May  1  image  (Fig.  1)  has 
Am<0”  1.  The  departures  of  Am  from  zero  are  probably  an 
artifact  of  our  preliminary  calibration  of  V2,  which  can  be 
biased  by  calibration  errors.  The  shape  of  the  <-/>cl(X)  curve 
and  the  size  of  the  jump  are  much  less  sensitive  to  calibration 
effects. 

4.  DISCUSSION 

4. 1  Broadband  Multichannel  Data 

The  images  shown  here,  with  a  spatial  resolution  of  ^  3 
mas  and  generated  from  five  to  eight  minutes  of  data  each, 
show  the  power  of  broadband  multichannel  optical  interfer¬ 
ometry.  Simultaneous  multichannel  observations  over  a 
broad  band  produce  data  at  a  wide  range  of  spatial  frequen¬ 
cies  even  with  only  three  baselines,  providing  a  dramatic 
improvement  in  (u,v)  coverage  over  that  of  a  monochro¬ 
matic  three-element  array  (e.g.,  Baldwin  et  al.  1996).  Optical 
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Fig.  4.  Images  of  Mizar  A  on  1996  May  2,  3,  4,  29,  and  June  1  and  4  (panels 
1-6  respectively).  Contour  levels  are  —  2%,  2%,  5%,  10%,  20%,  30%,  40%, 
60%,  and  80%  of  the  peak.  The  3  mas  (FWHM)  circular  Gaussian  restoring 
beam  is  shown  in  panel  1. 


interferometers  can  perform  well  over  a  wide  wavelength 
range,  while  radio  astronomical  receivers  have  been  limited 
for  historical  reasons  to  a  number  of  narrow  discrete  bands  at 
wavelengths  spaced  by  factors  of  roughly  two. 

Future  optical  interferometric  imaging  will  continue  to 
depend  on  large  total  bandwidths  to  obtain  good  ( u,v )  cov¬ 
erage.  For  a  number  of  practical  reasons,  such  as  the  fact  that 
low-noise  coherent  optical  amplification  is  not  possible,  the 
maximum  number  of  elements  in  an  optical  interferometer  in 
the  foreseeable  future  is  less  than  ten.  The  imaging  sub-array 
of  the  NPOI,  with  six  elements  when  complete,  will  have  15 
baselines,  as  compared  with  the  351  baselines  of  the  VLA. 
However,  each  of  the  NPOI  baselines  will  take  data  in  32 
channels  spanning  XX  450—850  nm,  thereby  expanding  its 
(u,v)  coverage  by  nearly  a  factor  of  two  and  increasing  the 
number  of  spatial  frequencies  to  480.  For  objects  in  which 
the  image  structure  is  not  independent  of  wavelength,  the 
analysis  must  generate  images  in  an  ( a,S,\ )  cube  rather 
than  in  (a,  8)  alone,  but  the  advantage  of  a  large  number  and 
range  of  spatial  frequencies  is  maintained. 


4.2  Vacuum  Delay  Lines 

The  vacuum  delay  fines  of  the  NPOI  are  crucial  for  ob¬ 
taining  data  over  such  a  wide  band,  particularly  for  extend¬ 
ing  the  band  to  the  blue.  The  greatest  part  by  far  of  the 
differences  in  path  length  from  a  star  to  the  interferometer 
apertures  is  in  vacuum,  above  the  Earth’s  atmosphere.  In  the 
interferometer,  these  vacuum  path  length  differences  can  be 
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Fig.  5.  Binary  separations  and  uncertainty  ellipses  as  obtained  from  the 
NPOI  data,  superposed  on  the  orbit  of  Hummel  et  al.  (1995)  from  Mark  III 
optical  interferometer  data.  The  separations  are  from  fits  to  the  images  in 
Figs.  3  and  4.  The  ellipses  are  one  fifth  of  the  synthesized  beam,  a  standard 
VLBI  indication  of  the  uncertainty.  The  measured  separations  are  connected 
by  lines  to  the  separations  predicted  from  the  Mark  III  orbit.  The  major  axis 
of  the  orbit  is  shown,  with  “T”  marking  the  periastron.  The  spacing  of  the 
dots  along  the  orbit  corresponds  to  P/360,  where  P  =  20? 54.  A  small  arrow 
indicates  the  sense  of  the  orbital  motion. 


compensated  in  vacuum  or  in  air.  If  they  are  compensated  in 
air,  additional  dispersive  optics  may  be  needed. 

When  the  vacuum  path  length  differences  are  compen¬ 
sated  with  vacuum  delay  lines,  the  only  remaining  differ¬ 
ences  in  air  path  lengths  between  beams  are  due  to  atmo¬ 
spheric  turbulence,  so  that  fringe  smearing  due  to  the 
variation  in  longitudinal  dispersion  across  a  spectral  channel 
or  across  the  observed  band  is  reduced  to  manageable  pro¬ 
portions. 

If  one  used  air  delay  lines  but  no  dispersion-correcting 
optics,  and  if  n(a)~  1  of  air  varied  only  linearly  with  the 
wave  number  cr,  one  could  eliminate  the  variation  of  fringe 
phase  with  wavelength  by  suitably  adjusting  the  delay  lines. 
But  because  of  the  higher-order  terms  in  n{cr)—  1,  which 
become  significant  in  the  blue  half  of  the  visual  band,  elimi¬ 
nating  the  variation  of  phase  with  wavelength  in  this  way  is 
not  possible.  As  an  example,  for  10  m  of  air  path  mismatch 
between  beams,  setting  the  delay  lines  to  minimize  the  varia¬ 
tion  of  fringe  phase  still  leaves  >170  radians  of  variation 
within  the  XX  450-850  nm  band  and  more  than  5  radians 
of  variation  within  most  of  the  32  individual  channels  of  the 
NPOI  spectrometers. 

For  relatively  short  baselines,  one  can  use  air  delay  lines 
and  compensate  for  the  atmospheric  dispersion  with  glass. 
With  the  appropriate  amount  of  a  typical  glass  such  as  BK7 
in  one  interferometer  beam  compensating  for  1 0  m  of  air  in 
the  other  beam,  the  residual  variation  of  fringe  phase  across 
the  XX  450—  850  nm  band  is  —0.3  radian.  But  for  baselines 
of  100  m  or  more,  and  the  corresponding  air  path  mis¬ 
matches,  the  residual  phase  variations  after  correction  are 
still  severe.  Even  the  variation  over  a  10  nm  spectral  channel 
is  as  large  as  0.25  wave,  large  enough  to  lead  to  greatly 
reduced  fringe  visibilities. 

The  large  total  bandwidth  of  the  NPOI,  and  therefore  the 
use  of  vacuum  delay  lines,  was  originally  driven  by  the  need 
for  milliarcsecond  astrometric  accuracy.  Using  a  large  band¬ 
width  enables  us  to  measure  the  variation  of  fringe  phase 
with  wavelength,  estimate  the  air  path  mismatch  from  this 
variation,  and  remove  the  atmospheric  contribution  to  the 
delay.  The  remaining  geometric  delay  is  the  fundamental  as- 
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trometric  datum  for  an  interferometer.  The  results  presented 
here  show  the  utility  of  large  bandwidths  for  imaging  as 
well. 

4.3  New  Analysis  Techniques 

In  radio  interferometry,  the  fundamental  data  are  the  com¬ 
plex  visibilities  Ve from  pairs  of  elements.  At  optical 
wavelengths,  however,  the  phase  of  the  complex  visibility  is 
completely  dominated  by  atmospheric  fluctuations,  which 
have  a  timescale  of  ~  10  ms  for  one-radian  fluctuations;  thus 
the  phase  carries  no  information  about  source  structure.  In 
addition,  the  obvious  visibility  amplitude  estimator 
V=(C/N)(X2+Y2)112  [where  C  is  the  normalization  factor 
in  Eq.  (1)]  is  biased,  since  (X2+Y2)  is  nonzero  (even  for 
zero  visibility)  due  to  detector  noise  (Tango  &  Twiss  1980), 
and  this  bias  is  a  function  of  V  itself. 

The  closure  phase,  however,  is  free  of  atmospheric  ef¬ 
fects,  and  the  V 2  estimator  of  Eq.  (1)  is  unbiased.  For  these 
reasons,  the  fundamental  data  for  optical  interferometric  im¬ 
aging  are  the  V 2  from  the  interferometer  elements  taken  in 
pairs  and  the  4>c\  for  the  elements  taken  in  triples.  In  cases 
where  the  visibility  on  only  one  of  the  three  baselines  of  a 
triple  is  small,  the  amplitude  of  the  triple  product,  V,  V2V3 , 
will  have  a  higher  signal-to-noise  ratio  than  V2  for  that  base¬ 
line. 

Our  results  demonstrate  the  need  for  analysis  techniques 
that  are  suited  to  V2  and  <f>d  data.  We  used  the  standard 
radio-interferometric  techniques  for  our  images,  as  did  Bald¬ 
win  et  al.  (1996)  for  their  Capella  images,  by  using  the 
square  root  of  V2  and  by  assigning  tpci  to  one  of  the  base¬ 
lines.  But  using  these  techniques  has  two  disadvantages.  One 
is  that  (V2)112  is  biased:  it  is  always  positive,  but  the  mea¬ 
sured  V 2  may  be  negative  due  to  noise  in  X  and  Y.  (This  bias 
does  not  appear  to  be  the  cause  of  the  non-zero  Am  in  our 
images,  however.)  The  other  is  that  we  cannot  adjust  the 
relative  weights  of  the  <pci  and  the  V2  data,  since  they  are 


treated  together  as  a  complex  visibility.  Adjusting  these  rela¬ 
tive  weights  is  desirable,  since  in  general  they  have  different 
signal-to-noise  ratios.  In  the  case  of  our  data,  the  shape  of  the 
jump  in  (f>cl  is  the  more  sensitive  indicator  of  Am. 

5.  CONCLUSION 

We  have  presented  the  first  images  from  the  Navy  Proto¬ 
type  Optical  Interferometer  (NPOI).  Multichannel  observa¬ 
tions  over  a  large  total  bandwidth  provided  broad  ( u,v )  cov¬ 
erage  even  though  we  used  only  three  baselines.  The  vacuum 
delay  lines  of  the  NPOI  make  it  possible  to  avoid  the  effects 
of  longitudinal  dispersion  over  a  large  band.  We  processed 
the  data  with  standard  radio-interferometry  techniques  to 
produce  our  images.  However,  the  fundamental  observables 
of  optical  interferometry,  V2  and  <f\\ ,  are  different  from  the 
complex  visibilities  Ve of  radio  interferometry.  Analysis 
techniques  that  are  more  appropriate  to  the  optical  interfer¬ 
ometry  case  are  clearly  needed. 

The  resolution  of  mas  in  our  images  of  the 
20'?54  spectroscopic  binary  f 1  UMa  (Mizar  A)  is  the  highest 
ever  attained  in  an  image  at  visual  wavelengths,  and  exceeds 
by  a  factor  —2  —  3  the  best  resolution  theoretically  obtain¬ 
able  by  the  largest  current  monolithic  optical  telescopes.  The 
separations  and  position  angles  of  the  components  of  Mizar 
A  in  our  images  show  good  agreement  with  the  9.6  mas  orbit 
measured  with  the  Mark  HI  optical  interferometer  on  Mt. 
Wilson  (Hummel  et  al.  1995),  but  show  small  systematic 
differences  with  the  predictions  from  the  Mark  III  orbit.  The 
differences  may  be  due  to  a  gap  in  the  orbital  phase  coverage 
in  the  Mark  III  data.  The  closure  phase  data  are  very  sensi¬ 
tive  to  Am;  our  preliminary  estimate  is  that 
Am~0”01±0“01. 

This  research  was  funded  by  the  Office  of  Naval  Research 
and  the  Oceanographer  of  the  Navy. 
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